SINGULAR ASYMPTOTIC EXPANSIONS FOR DIRICHLET 
EIGENVALUES AND EIGENFUNCTIONS OF THE LAPLACIAN 
ON THIN PLANAR DOMAINS 

DENIS BORISOV AND PEDRO FREITAS 

Abstract. We consider the Laplace operator with Dirichlet boundary con- 
ditions on a planar domain and study the effect that performing a scaling in 
one direction has on the spectrum. We derive the asymptotic expansion for 
the eigenvalues and corresponding eigenfunctions as a function of the scal- 
ing parameter around zero. This method allows us, for instance, to obtain 
an approximation for the first Dirichlet eigenvalue for a large class of planar 
domains, under very mild assumptions. 



1. Introduction 

The study of the spectrum of the Laplace operator on thin domains has received 
much attention in the mathematical literature over the last few years. Apart from 
the connection to certain physical systems such as quantum waveguides, this limit 
situation is also of interest as it may provide insight into certain questions re- 
lated to the spectrum of the Laplacian - see the recent paper by Friedlander and 
Solomyak |FS| for some references on both counts. 

The purpose of the present paper is to study this spectrum in the singular limit 
around a line segment in the plane. More precisely, given a planar domain we 
consider its scaling in one direction, so that in the limit we have a line segment 
orthogonal to this direction. In particular, and under mild smooth assumptions on 
the domain, we derive the series expansion for the first Dirichlet eigenvalue thus 
showing that the coefficients in this series have a simple explicit expression in terms 
of the functions defining the domain - see Theorem [T] below. 

Due to the notorious lack of explicit expressions for Dirichlet eigenvalues of 
planar domains, this seems to be a possible path towards obtaining information 
about such eigenvalues, and indeed it was one of the motivations behind our work. 
As an example, consider the case of the family of ellipses centred at (1/2,0) and 
with axes 1/2 and e/2. In this case our results yield that 

(1.1) Ai(£) = ^ + — +3+ + j)e + 0(e2), as e ^ 0, 

thus complementing the asymptotic expansion for the first eigenvalue of ellipses 
around the disk derived by Joseph in [j]. Further terms in the expansion may 
be obtained by means of Lemma 13.21 The graphs of the four-term asymptotic 
expansions for two other examples, namely the lemniscate and the bean curve are 
shown in Figure[TJ For a discussion of these and other examples, see the last section. 
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Figure 1. Graphs of the four-term asymptotic expansions ob- 
tained for the first eigenvalue of the lemniscate and the bean curve 
in Examples [2] and [3] presented in Section [5] The points are nu- 
merical approximations to the eigenvalue with an error not greater 
than 10-5 . 

The advantages of our method are twofold. On the one hand, it does not require 
any a priori knowledge neither of the eigenvalues nor of the eigenfunctions of a 
particular domain, as the perturbation is always made around the singular case of 
the line segment. In the above example, for instance, Joseph's method requires 
knowledge of the first eigenvalue of the disk. This limits applications of the same 
type of argument to general examples. Furthermore, it is not evident that even in 
cases where the eigenvalue of the unperturbed domain is known one can determine 
explicitly the value of the coefficients of the perturbation around this domain in 
terms of known constants. This is illustrated by the example in Section 5 in [j], 
where one of the coefficients in the expansion gives rise to a double series. On the 
other hand, and by its own nature, our method is particularly suited to dealing 
with long and thin domains where numerical methods will tend to have problems. 

The fact that we consider a singular perturbation does pose, however, some dif- 
ficulties. In order to understand these, let's begin by observing that the behaviour 
of the asymptotic series obtained for the eigenvalues depends on what happens at 
the point of the line segment where the vertical width is maximal, say x. More 
precisely, the coefficients appearing in the expansion will depend on the value that 
the functions defining the domain and their derivatives take at x, but not on what 
happens away from this point. This is not surprising, since it is in this region that 
the first eigenfunction will tend to concentrate as the parameter goes to zero. In 
order to derive the asymptotics, we assumed that there existed only one such point 
and that the domain was C°° smooth in a small neighbourhood of x. While the for- 
mer condition may be somewhat relaxed, and a finite number of points of maximal 
height could, in principle, be considered, if the latter condition is not satisfied then 
the powers appearing in the asymptotic series will actually be different. 

A first point is thus that the type of expansion obtained will depend on the local 
behaviour at x. To illustrate what may happen, consider the following example 
studied by the second author using a different approach [F] . Let Tp be an obtuse 
isosceles triangle where the largest side is assumed to have fixed unit length, and 
where the equal angle P approaches zero. In this case we have 




2 4-2^/Vi^^/3 



+ 0(/3-2/3) as/3->0. 
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where a[ « —1.01879 is the first negative zero of the first derivative of the Airy 
function of the first kind. We thus see that having the maximal width at a corner 
point introduces fractional powers in the expansion. 

On the other hand, it can be seen from our results that what happens away 
from X does not influence the coefficients appearing in the expansion in Theorem [1] 
Thus in the example given above for ellipses, changing the domain outside a band 
containing the mid-region will yield precisely the same coefficients although the 
expression for the first eigenvalue as a function of e will certainly be different. This 
means that the series described in Theorem [T] cannot, in general, converge to the 
desired eigenvalue, and that there should exist a tail term which goes to zero faster 
than any power of e, which our analysis does not allow us to recover. An interesting 
question is thus to determine the nature of such a term and when can we ensure 
that the series expansion derived here is actually convergent to the corresponding 
eigenvalue. 

Note that although we address the same problem as in |FS| , both our approach 
and results arc quite different from those in that paper. In particular, we obtain 
the full asymptotic expansion for the eigenvalues and eigenfunctions. Furthermore, 
the coefficients in these expansions are obtained as a simple explicit function of the 
value of the functions involved and their derivatives at the point where the function 
H takes its maximum. Another difference is that we consider a two-parameter set 
of eigenvalues, while the method used in |FS| leads to a one-parameter set. On 
the other hand, our approach required the function H to be smooth in a neigh- 
bourhood of this point, while the results in [FSj are more general in this respect, 
allowing for the existence of corners, for instance. We hope to be able to consider 
this situation in a forthcoming paper. Furthermore, in |FSj the authors actually 
prove convergence in the norm of the resolvents, which is a stronger property than 
convergence of the eigenvalues and eigenfunctions. 

The plan of the paper is as follows. In the next section we state the results in 
the paper. In Section[3]wc prove the general form of the asymptotic expansions for 
eigenvalues and eigenfunctions. The proof is split into two parts. First we construct 
the asymptotics expansions formally. The main idea here is to use the ansatz of 
boundary layer type which localizes in a vicinity of the point x mentioned above. 
After formal constructing, the formal expansions are justified, i.e., the estimates 
for the error term are obtained. Here the main tool is Vishik-Lyusternik's lemma. 
Section m contains the study of the expansion for the first eigenvalue. This consists 
in identifying the eigenvalue which corresponds to a positive cigcnfunction, and 
then obtaining the explicit form of the terms in the expansion of the eigenvalue. In 
the last section we consider some applications of our results to specific domains. 



2. Statement of results 

Let h± = h±{xi) G C[0, 1] be arbitrary functions and write Hixi) := /i+(a;i) -I- 
h-{xi) for xi in [0, 1]. We shall consider the thin domain defined by 

Vl^ = {x : xi ^ (0, 1), — £/i„(xi) < X2 < eh+[xi)}, 

for which we assume that the function h attains its global maximum at a single 
point X e (0,1), and that H{x{) > for xi € (0,1). Note that the cases where 
either H{0) = or H{1) = are not excluded. We also assume that the functions 
h± are infinitely differentiable in a small neighbourhood of x, say {x — d,x). For 
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the sake of brevity, in what foUows we shall write Hq := H{x) and denote by Hi 
the i-th derivative of ft. at x. In the same way we denote the derivatives of h- by 
hi. We shall assume that there exists fc ^ 1 such that 

(2.1) H, = 0, i-l,...,2fc-l, H2k<0. 

Our aim is to study the asymptotic behaviour of the eigenvalues and the eigen- 
functions of the Dirichlet Laplacian — Aq^ in il^. Let x = x{t) G C°°(K) be a 
non-negative cut-off function which equals one when |^i — a;| < 5/3 and vanishes 
for \t-x\> S/2. We denote := n {x : \t - x\ < 6}. 

The main results of the paper are the following two theorems. 

Theorem 1. Under the above conditions, there exist eigenvalues A„.„i(e), n,m € 
N, of the operator — whose asymptotic expansions as e goes to zero read as 

oo 

(2.2) A„.„(£)=rt^'")+e-2^,7M"■"^ V-=e", " fcTT' 

i=2k 

where 

2 2 

(n.m) 7T 71 (n.m) * (n.ni) ^ 

''0 = tt2 ' '^2k ~ ^^n,m, — U, 

("^™) ^^"^^1 27r^n^g2fc+l ^^2fc+l.T,(",'") ^ \ 

(2.3) C,,^2 = - (2fc+l)!g3 (gl *1 ^ *r,„OL.(M) 



2^^n^g2fc+2 ||^fc+i , . 3^^n^g|fe ^2fc^ „ 

(2fc + 2)!H3"^i ((2fc)!)'i/4ll^i <P".™llL.(E)- 

_ffere A„.„i bkc? $ri.m are the eigenvalues and the associated orthonormalized in 
L2{R) eigenfunctions of the operator 

- de, {2k)\Hi 

in L2{M.), Sik is the Kronecker delta, and xj/^"'™'' g C°°(]R) n L2(K) is the exponen- 
tially decaying solution to 

which is orthogonal to <&n,m in ^2(1^). 

Given an eigenvalue Xn^mi^), let Al*m(e) be the eigenvalues o/— A^ having the 
same asymptotic expansions 112. S^) as A„.m(e), and Xn)mi£) '■— A„.m(e). Then there 
exists a linear combination 4'n,m{x,e) of the eigenfunctions associated with \n]m(£), 
whose asymptotic expansion reads as follows 



00 



(2.4) ^„,„.(a;,£)=x(xi)V77^^i"''"HO, 



i=0 



in the W2{^e)-'^orm and W2{^i^^)-norm, where 
_ xi - X _ X2 + sh^{xi) 

~ eH{xi) ' 

= 1>n,™(6)sin^<2, V^"'™^(0 = *i"'™^(6)sin7r<2. 
The remaining coefficients of the series 112. ^2.4^ are determined by Lemma 
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The second result gives explicit expressions for the first four nonvanishing terms 
in the asymptotic expansion of the first eigenvalue in terms of the functions h± and 
their derivatives at x in the case where H2 is negative. 

Theorem 2. For any N ^ 1 there exists Eq ~ £o{N) > such that for e < Eq 
the first N eigenvalues are Ai_m(£), m = 1,. . . ,N. These eigenvalues are simple, 
and the asymptotic expansions of the associated eigenf unctions are given by \2.4-^ , 
where tjji.m stands for the corresponding eigenfunction. In particular, if k ^ 1, the 
lowest eigenvalue will have the expansion 



where 



Ai,i(e) 


(14; 

Co 


-+ 2 


-+ci^^^^ 








(1.1) 
Co = 


-Wo 


(1.1) 
C2 = 


7r( - H2) 

Hq 


(14) 


2 1 2 
TT hi 


9i/2 


UHi 


Hi 


16i?o 





f 4^'^^£ + C'(£2), as£^0, 
1/2 



(2.5) 



c, 



(1,1) ^ hI'^ U^hl TT^{hiH3+Hl)h2 83g| 
' -^-if2)'^' V^^" ^^0^2 256^2 

~3MHf ~ 691277| 384H^ ~ USH^ ^ 576iJ| 
7i?| i/e TT^Hf TT^ffa/i? TT^hshA 



768iJ| 192^2 6iJ^ 2i/3 2i?2 |- 

In the case of higher eigenvalues, which of the eigenvalues A„,,„(£) corresponds 
to a given eigenvalue will depend on the value of s. Furthermore, generally 
speaking, the eigenvalues given in Theorem [T] do not exhaust the whole spectrum. 
Indeed, assume that the function H has a local maximum at a point x ^ x. Then 
one can reproduce the proof of Theorem[T]and show that there exists an additional 
infinite two-parametric set of the eigenvalues associated with x; the corresponding 
eigenf unctions are localized in a vicinity of x. 

Let 7 := Ai,2 — Ai,i denote the length of the first spectral gap. The above result 
allows us to obtain the first term in the asymtotic expansion for 7. 

Corollary 2.1. Under the above hypothesis, the quantity £^{£) remains bounded 
as £ goes to zero. If k = 1, then 

7(e) = ^^^^§^+0(1), ase^O. 

Hq £ 

We remark that, in general, the spectral gap is unbounded when either the di- 
ameter or the area are kept constant, the simplest example being that of a circular 
sector where the opening angle approaches zero. From the results for obtuse isosce- 
les triangles in [F] we know this to be the case also for one-parameter families of 
domains of the type considered here. What this result shows is that under the 
conditions above the gap will remain bounded as e goes to zero if we keep the area 
fixed. To see this, it is sufficient to note that \^e\ = £|i^i| and thus |r2e|7(£) is also 
bounded. In particular, this shows that regularity at the point of maximum height 
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plays an important role in bounding the gap. Note that this boundedncss is not 
uniform on the domain. Also, if we fix the diameter instead of the area then the 
gap will still go to infinity as e goes to zero. For a sharp upper bound for the gap 
and a numerical study of the same quantity, see |FK| and [AF| , respectively. 



3. Proof of Theorem [T] 



Let Ag and ipe be an eigenvalue and an associated eigenfunction of — , re- 
spectively. We construct the asymptotics for and ips as the series 

oo oo 
(3.1) Ae = e^^/ie, Me = Co + ^ T]''ci, Ve(C) = ^ V'AiO, 

i=2k i=0 

where Ci and ipi are the coefficients and functions to be determined. 

In what follows we will show that the function ip^ is exponentially small (with 
respect to e) outside f^f. This is why we are interested only in determining ipi on f^f. 
After passing to the variables ^, the domain 17^ becomes : |^i| < Si]'^ , < ^2 < 
1}. As 77 0, the latter domain "tends" to the strip 11 := : < ^2 < !}• This is 
why we will consider the functions tpi as defined on 11. The mentioned exponential 
decaying of the eigenfunction is implied by the fact that all the coefficients ijji decay 
exponentially as |^i| — s- ±cx), in other words, we postulate the latter for ■i/'i. 

Having the made assumptions in mind, we rewrite the eigenvalue equation for 
■0£ considered in fi^ in the variables ^, 

V'e = on 911, 

where Kij = (^177, ^2, ^» = ^i(6?7, 6, '7), and 



2fc I. _ ^^^^^h'_ix + Zi) - Z2H'{x + Zi) 



H{x + zi) 

, , l + V^''+^h^X + Z,)^Z2H'{x + Z^)f 
K22(Z,1j} = -J— ^ , 

(3.3) h^(x + zi) 

2k-\-2 

K2{z,v) = -rrh ^(hlix + zi)H{x + zi) ~ 2h'_ix + zi)H'ix + zi) 

n^lx + zt ] \ 



+ 2z2{H'{x + zi)y - Z2H"{x + zi)H{x + zi)) . 



Now we expand the coefficients K12, K22, and K2 into the Taylor series with respect 
to 77, taking into account (|2.ip . 



if 12(6^,6,^)= E '?'(^i2^(a)+6Qi2'(6)) 

(3.4) K22{^iV,^2,v) = (^22^(6) +6^22(6) +^2^2(6) 

00 



i=2fc+2 
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where PI2, Pi-n Ph Qi2j Q221 Ql certain polynomials with respect to ^1, and, 
in particular, 

p(2fe+l) _ hi p(2fe+2) _ ^2 ^ 

^12 ~ Tj ' ^12 ~ rr '^1' 

^0 -no 

^(2fe+l) _ „ ^(2fc+2) _ H2 



(3.5) 



p(2A;) _ 2H2kS,l'' p(2fc+l) _ _2il2AH-l|r^ 

22 - (2fc)!773' 22 - (2fc + i)!i/3' 

p(2fc+2) _ '^H2k+2iT^'^ I c ^H^k^f' hi 

- (2fc + 2)!i?„3 +'^^-((2fc),)2^4 + i?2 

^(2fc) _ ^(2fc+l) _ ^(2fc+2) _ „ p(2fe) _ p(2fe+l) _ p(2fe+2) _ „ 
^^22 ~ ^^22 ~ ^22 ~ ^' ^22 ~' ^22 ~ ^22 ~ ^' 

p(2A:+2) _ h2_ n(2*:+2) _ -^2 
2 — rr ' ^2 ^ tt ' 

We substitute (|3.ip and p.4p into (|3.2p and evaluate the coefficients of the same 
powers of 77 taking into account p.Sp . It leads us to the system of the boundary 
value problems, 

(3-6) TTl^^ + co^^^O in n, 1^, = on an, i = 0,...,2fc-l, 

1 d^^P2k n^n^ , 32.^0 2H2k^d^i^o , , • tt 
V'2fe = ^^72 7TT7T777T ^t72~ + C2fci/'0 m 11, 



(3.7) 8^2 Hi "^^^ {2k)\Hl 9^1 

■02^ = on an, 
1 a2V'. , , , dH\ 



CQlpi =Ci1po + i^'n \- C2ki^i~2k 



Hi de2 -^'^^ ■ dii 

(3.8) 2i72feC?' ^H^^-2k 



in n. 



(2fc)!i/3 d^2 
■04=0 on an, z^2fc + l, 
where the functions ipQ, tpi are assumed to decay exponentially as ^1 — > ±00, and 



2 



F,=2 ^ (p!^'^(a)+6Ql^'(e 

i=2fc+l 



j=2fe+l ^2 



E (^i2'(ei) + 60^^'(a) + e2'i?i'(ei)' 

-1 



E (^i^'^(ei)+6Q^^'^(a))^+ E 



9^ 

j=2fc+2 ^ j=2A:+l 



Problems (|3.6p can be solved by separation of variables. It gives the formulas 
for tpi^ i = 0, . . . ,2k ~ 1, and cg, 

2 2 

(3.9) 7/',(0 = *»(6)sinri7re2, cq = — 5-, n € N, 

where i = 0, . . . , 2fc — 1, and the functions \E'i are to be determined. We can consider 
the problem p.7p as posed on the interval (0, 1) and depending on the parameter ^i. 
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Hence, this problem is solvable, if the right hand side is orthogonal to sin7T,7r^2 in 
^2(0, 1) for each ^1, where the inner product is taken with respect to ^2- Evaluating 
this inner product and taking into account p.9p . we arrive at the equation 

(3.10) (-^- i2kyHt J""^ " 

The exponential decaying of -00 a-s |^i| — > 00 is possible, if the same is true for 
^0- Hence, ^"0 is an cigenfunction of G^, and therefore. C2k is the corresponding 
eigenvalue. Thus, = ^n,m, and the formula (|2.3|) for C2k is valid. 

Lemma 3.1. The spectrum of Gn consists of infinitely many simple positive iso- 
lated eigenvalues An^m- The associated eigenf unctions are infinitely differentiable, 
and decay exponentially at infinity, namely, 

'i>„,m(ei) = C^±iai^"^e-^l«^l'(l + o(l)), if k^l, 

i>n.™(ei) = c'±iar*e-^i«^i'^\i+o(i)), if k>i. 

The last formulas can be differentiated with respect to ^1 . 

The equation (G„ — Kn,m)u = f , f (z L2{R), is solvable if and only if 

(3.11) (/,$n.m)L.(M) -0. 

In this case, there is a unique solution u orthogonal to in L2{S.). If f (z G°°(]R) 

is an exponentially decaying function satisfying 

(3.12) / = 0(161'' exp-^l«^l'^'+'), 161 -ex., 

and this identity can be differentiated with respect to 6; then the solution u is 
infinitely differentiable and decays exponentially, 

(3.13) ^i = 0(|6|'^exp-^l«^l'"^'), |6I - ex., 

where [3 is a some number. This identity can be differentiated with respect to 6. 

Proof. The statement on the eigenfunctions follows from |BSh| Ch. I, Sec. 2.3, Th. 
3.1, Sec. 2.4, Th. 4.6]. The solvability condition (j3.1ip is due to self-adjointness 
of Gn- Theorem 4.6 in |BSh[ Ch. E, Sec. 2.4] gives also the asymptotic behaviour 
of the fundamental system of the equation (j3.10p . Using these formulas and repre- 
senting the solution u via Green function, one can easily prove (|3.13p . □ 

Taking into account (|3.9p . p.lOp . we can rewrite (|3.7p as 



1 a^V': 



" V'2A: = in H, Tp2k = on 9H. 



Hence, the formula (|3.9p is valid for ip2k as well. In what follows, and for the sake 
of brevity we will write simply A and $ instead of A„.,„ and '^n,m- 

Lemma 3.2. The problems iS. 6\) . 7| ), iS. 8]) are solvable, and their solutions read 
as follows: 

(3.14) V'^(0-v^«(c) + *.(a)sin7r<2, v^«(e) - ^ (^i)<^S fe), 

3 

where the sum with respect to j is finite, e G°°(M) are exponentially decaying 
functions satisfying ![3l3\) . cjyf] G G°°[0, tt] are orthogonal to sin7rn^2 in -^2(0, tt) 
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for each G M and vanish on 911. As z = 0, . . . , 2k, the functions are identically 
zero, while for i > 2k they solve the boundary value problems 



1 dH', 7r2n2~ d^^J,-2k 2iJ2fc^f'^ aV»-2fe 

-Ipi = 



(3.15) 



Hi de. Hi de, i2ky.Hl 

+ ATpi^2k + ^ CjV'i-j + Fi - fi sin7rn^2 in n, 

j=2k+l 

tpi =0 on dU, 



F,=2 ^ (p/^V^Qis' 



j=2fe+l 

32 



(3.16) + E (Pi^V6Q^^'+?l42;^ 



96 ' 



i=2fc+2 

(3.17) /i = 2(FJ,sin7rr^6)L2(o,7r)■ 
T/le functions £ C°° (R) satisfy i3.13\) . are orthogonal to $ m ^2 (M) and solve 
the equations 

i-l 

(3.18) (G„ - A)'i',_2fc = /. + E 

j=2fc+l 

The numbers Ci, i ^ 2fc + 1, are given by the formulas 

(3.19) c, = -(/,, $)l,(r). 

Proof. Wc prove the lemma by induction. The statement of the lemma for i = 
follows from (|3.9|) . This identity also implies the formulas (|3.14p for « = 1, . . . , 2fc, 
where i/'j = Oj and ^'i are some functions to be determined. 

Assume that the statement of the lemma holds true for i < p. p ^ 2fc + 1. Then, 
in view of p.l4p . the right hand side in the equation in (|3.8p can be rewritten as 

\ 

Cp'^ ~ {Gn - A)*p-2fc + E ^j'^P-J sin7r<2 

(3.20) ^ ' 

+ -W " miHi aci + ^^^^ + 

The solvability condition of p.8p is the orthogonality of this right hand side to 
siuTrn^ in -L2(0,7r) for each 6 € We write this condition, taking into account 
the orthogonality of ipj, j < p, to siuTrn^ in L2(0,7r), and the relation 



-2k 



„ ,sin7m6 = -TT^n^ ( V'p-2/c,sin7m6 , 

/ -L2{0,7r) 



0. 



This procedure leads us to (|3.18p . By Lemma [01 the solvability condition of (|3.18p 
is exactly the formula p.l9p . since the functions 5'i_2fc, * < P, are orthogonal to <I> 
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in L2(M). We choose the solution of this equation to be orthogonal to $ in L2(R) 
and note that by the formula (|3.14p for ^E';, i < p, the function fi satisfies p.l2|) . 
Hence, by Lemma [3.11 the function ^'p_2fe satisfies (|3.13p . The formulas p.l4p 
and (|3.20p yield that the right hand side of the equation in (|3.8p with i = p can be 
represented as a finite sum ^ f^]{ii)f'^] (^2), where /^^j € C°°[0, 1] are orthogonal 

to sin7r<2 in -^2(0,77), while the functions /^^j e C°°(R) satisfy ((37T2)l . This fact 
implies the formula (j3.14p for i = p. □ 

Let us prove the remaining formulas in (|2.3p . It follows from (|3.16p . p.Sp that 

(3.21) ^^2..+! = ^^cos,r<2 + (2fc+i),g3 -&^in7r<2. 

Therefore, 

9_2„2r7- f2fe+l 

(^■2^) = (2fe+i;!l3 ^- 

The eigenfunction $ is either odd or even with respect to ^1 , due to the evenness of 
the potential in the operator G„. Hence, the function <I>^ is even, and this is why 
(/2fc+i,$)L2(K) = 0, proving the formula for C2fc+i. 
Employing (|3.5p , by direct calculation we check that 

/Q O-JA f ^{2k+2) 2 2n{2fe+l),T, 2 2n(2fe+2)* l^(2/£+2)^ 

(3.23) /2fe+2 = -Qi2 ^^-'^"-^22 *i-7rnP22 ^-o'32 ^■ 
By integrating by parts we obtain 

Now the formula (|2.3p for C2fc+2 follows from the last two identities and (|3.5p . 

We proceed to the justification of the expansions p.ip . For any > 7{k + 1) 
we denote 

N N 
1=0 i=2k 

Lemma 3.3. The pair ■0s, ^e.N satisfies the boundary value problem 

(3.24) - ^4'e,N ^ K,N'4's,N + g^N in n^, i>e,N^O on dfle, 
and 

N o 

\\9e,N\\L2{i\) = 0(e^~^). 

This lemma follows from the definition of the problems p.6p . (|3.8p . Lemma [321 
and the exponential decaying of ^/j^ as |^i | — > oo. 
We can rewrite the problem p.24p as 

(3.25) 1pe,N = iXe,N + l)Ae1pE,N + Aege,N, 

where indicates the inverse of — A^^ + 1. It is clear that the operator A^ is self- 
adjoint, has a compact resolvent, and satisfies the estimate ^ 1, uniformly in 
£. We rewrite (|3.25p as 

^E,N + t « + i 
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It follows from Lemma [3.31 and the estimate WA^W ^ 1 that 

1 

Aege,N 



where Cn are some constants. We apply Lemma 1.1 in |OSY| Ch. n, Sec. 1.1], 

where we let a = Cns'^''' (see inequality (1.1) in [OSYl Ch. 11, Sec. 1.1]). This 

lemma yields that given N for each e there exists an eigenvalue TAr(e) of such 
that 

\TN{e)-{l + Xe,Nr'\=0{e^-^). 
Hence, there exists an eigenvalue Ajv(£) of — A^^ such that 

(3.26) |Ajv(e) - A,^Ar| Cns^-\ 

where Cn are some constants. 

Let En be a monotone sequence such that Cne''^ ^ Cn-i as e ^ Sn- Given n 
and m corresponding to the series (|3.1[) (see (|3.9p . (|3.10|) ). we chose the eigenvalue 
A„,m(e), letting A„,m(e) A^vle) as e e [eN,eN+i)- Employing p.26p . one can 
check easily that the eigenvalue Xn,m satisfies (|2.2p . 

For £ S [eN,£N+i), we employ Lemma 1.1 in |0SY1 Ch. IH, Sec. 1.1] once again, 
and we let in this lemma a ~ Cate"'^"^, d ~ 2Cn£''^~^ ■ It implies the existence 
of a linear combination tpn,m described in the statement of the theorem such that 

(3.27) - ^e,w||L.(o.) = OieT^-^) 
for each N. 

We denote (i := Xi, C2 := X2e-\ 1^ {C : Ci e (0, 1), -/i_(Ci) < C2 < /i+(Ci)}, 

(3.28) i'NiO ip7i.rn{x,e) - V'e,w(a;). 
Lemma 3.4. The relations 

l|V^/'Ar|U,(o,) =0(£^<'=+" =), IIV^JvIIm/KOVS) =0(£^(''+" =) 

/loM true. 

Proof. Employing Lemma 13.31 and (|3.27[) , (|3.26[) , by integrating by parts one can 
easily check that 

^ N 7 

(3.29) liVV'jv||L,(o.) =0(e^"^). 

It follows from Lemma [5751 the equation for ipe,N and (|3.27p . (|3.26p . (|3.29p that the 
function (/)((^) := xixi)ipn,mix, e) — ipe^Nix) satisfies the boundary value problem 

(3.30) - (^e^^ + ^^^^g in , = on dn' , 
where 

N 9 N Q 

(3.31) \\9\\L,in^)=Oie^-^), ||V0|U,(o^) = 0(eWT)-3). 

Integrating by parts in the same way as in the proof of Lemma 7.1 in [LUl Ch. 3, 
Sec. 7], one can check that 
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Here := dn^ \ {C : Ci = S} \ {( ■ Ci = S}, = i^is), v = {vi,V2) is the outward 
normal to and s is the arc length of . Employing the identity 



we continue the calculations 



pa 



dCi 



V2 



if = on ^^ 
OS 



d 



ds = i/i — — T^ds= / i^iTr^^2T;-ds 

oiy OS ou 



pa 



d fdcj) 



^1^2 ( — ) ds + - / ( ^ ) ds = - / (l^iJ/^ 



9s 9^2 

2 



90 9 90 



ov 



pa 



where ' denotes the derivative with respect to s. The obtained formula implies that 



90 
9a 



92 



92 



"^'9(1 '^9Ci9C2 



ds 



Ce- 



where the constant C is independent of e. Hence, by l|3.32p 
'9^0 920\ 9^0 2 



dQl ' 9C| 



L2(0'') 

9^0 2 



9Ci9C2 



■^'ll'/'CcllL(n^) 



Ce- 



9Cf 



i2(n<5) 



92 



9CI 



L2(0^) 



+ 



92 



Employing this estimate, we obtain 

,|2 ^ ^''^ ' 



9C? 



2e' 



920 920 
9CF'9<f 



9Ci9C2 



L2(0'5) 



H'2i(n'5)' 



2 

i2(n'S)' 

920 2 

L2(n«) 



9CI 



□ 



^e4(i_2£)||0cclli,(oa)-Ce-i||0||2^(,,a). 
Combining this estimate with (j3.3ip . we complete the proof. 
The proven lemma implies the asymptotics (|2.4p . 

4. Proof of Theorem [2] 

We begin by showing that for sufficiently small e the first N eigenvalues are 
^i,m(e), m = 1, . . . , iV. If /i_ = 0, this statement follows from Theorem 1 and the 
arguments of Sec. 6.1 in [FSj . If /i_ is not identically zero, we cannot apply directly 
the above mentioned results from [FS| . but it is possible to extend their proof to 
the case /i_ ^ with minor changes. Below we list the required changes and refer 
to |FSj for the detailed proof. 

The first change is that in our case by the function h in |FSj one should mean H . 
Suppose for a while that H is strictly positive. The space is defined as consisting 
of the functions 



■ sm 



■n[x2 + eh^{xi)) 



eH{xi) eH{xi) 

The function v{x) in the potential defined in |FS[ Eq. (1.5)] should be intro- 
duced as 

n^{H'ixi)f {27rh'_ixi)-H'{xi)Y 



v{xi) 



•m^{xi) 



4:H^{Xl) 
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The number M, as in |FS| . should be Hq- The definition of the operators Qe and 
Qe remain the same as in [FS| . and the operator H is our operator G„. Under such 
changes Theorems 1.2, 1.3 in jFSj remain true. Throughout the proofs of these 
theorems the function sin should be replaced by sin ^^^^^hi^^^)^^^^ and in all 
the integrals w.r.t. X2 the limits of the integrations are — £ft._(xi), eh^(xi). The 
other arguments in the proofs of Theorems 1.2, 1.3 remain unchanged. Thus, these 
theorems are valid for the case /i_ ^ under the additional assumption H > 0. 

Employing the proven Theorems 1.2, 1.3 for the case /i_ ^ and proceeding as 
in [FS| Sec. 6.1] one can check that in the case H ^ the eigenvalues Ai,,„(e), 
TO = 1,...,7V, are the first eigenvalues of — A^^ for e small enough. We also 
conclude that these eigenvalues are simple. Hence, for each m = 1,...,N the 
eigenvalue Ai_m(e) is the unique eigenvalue of — having the asymptotics (|2.2[) for 
n = 1 and given m. Hence, the linear combinations V'l.rn introduced in Theorem [T] 
are the eigcnf unctions associated with the eigenvalues Xi.m- 

It remains to prove the formulas (|2.5p . In the case considered n = m = fc = 1. 
The operator Gi is the harmonic oscillator, and its eigenvalues and eigenfunctions 
are known explicitly. Namely, 



(4.1) 



Ai,i=0, $i,i(6) 



3/2 



The first identity proves the formula for Cj^'^-*. The formula for Cg^'^'' follows from 



It is easy to check by direct calculation that 



(14) 



1 TT^H, 



(4.2) 



l«$i 



18 



a(ee? + 3)<i>i 



3 

4^' 



11 TT^Fg 

48 614 ij3 ■ 



Now we substitute these identities into formula (|2.3p for C2k+2 and arrive at formula 
(1^ for 4^'^^ 
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To proceed further, we need the formulas for the coefBcients of the series p.4|) 
up to the order rf . They read as foUows, 
(4.3) 

(5) _ h'jHq - hiH2 ^2 p(6) _ KHq - hiH^ - 3/12^2 ^3 



.(5) _ J^yi n(6) _ HjHo - 3g| 



(5) _ 2/li^ -^2^3 ^5 _ -^5 ^5 

(6) _ 45H4H2H0 - mOHl - HeH^ + 30g|go , hshiHp - KjH^ + /ijgo ^2 

(5) _ 2/11^26 ^(6) _ feiiJs + 2/^2 g2 ^2 n(5) _ n _ ^1^? 

^22 ^ ' ^22 ~" ^2 "^l' -^22 ~ -^22 ~ ' 

(5) _ h3Ho--2hiH2 (6) _ ^4^0 - 5/12^2 - Sfeiffa ^2 

2 - i/2 4l' ^2 - 2i?o2 

^(5) _ _H3^ M6) _ 5g| - ^4^0 ^2 

^^2 — TF^i' ^2 — ;7772 
It foUows from (jXTSl) . ([S^, (jX^ that 

(4.4) V^(^'^'(0 = ii?o/ii(l- 26)sin^6- 

In view of (|3.23p . (|3.5p we can also find the function ^'2^'"'^'' explicitly, 
(4.5) 

_ f -'Hi ^2^4 



*9 = T^Ci + {9H4,H2Ho - llHiHo - 81H^ 



^ 288HoH^ ^1 ^ 128^2^ ^ M56H^e ^ USHod J ^' 

Now we use the formulas ([21]) for cf i 4, and (HH), ([Sll), ([3T6l) . ([STT]) . ([4T|1 . 

g21), gSl), dMl), (1131) and obtain 

(4.6) 

( - iS - 25i§|^ ( - ^^^^^"^^ + ^^^^2 + UHlH. 
-HI 



4320i/|ifo 



lOSmHiHie (5175g|go - 576712/1? /j22 + 81H^H2H^ - imulHl 



I2H2HI 



2TT'^H3hi - 2ATT'^H2hih2 + 'SH3H2H0 + Un^Hihi ) I $1,1 



The function /s is odd with respect to ^1, and therefore, in view of (|3.19p and 

^5 



evenness of $1,1, the formula (|2.5p is valid for Cg^'^''. Employing the obtained 
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(1,1). 



identity for /s, we can solve explicitly the equation ()3.18p for 'I' 
^(14) ^ _ ^^Hl 9 ^'^1-^3 , . . _ 9H4H2H0 



2f5 



64800i/(JiJ|6i V 

t3 



311040gog| (4455^4^3^2^0 - 2515H|i?o - SSlbH^H^ ~ mSHlH^Ho 



103680i72i7|6i 



1855i?|i7o + 3915i?4i?3i?2i^o + UU^HsH^Hq 



- 1728H5H^H^ - 1036807r^i/|/ii/i2 + 518407r^i/|/ii j j $1,1 
We substitute the relation obtained, the formulas (|2.5p for c-^'^'', i ^ 5, (|3.5p . p.9p . 



(IXTBl) . (IXT71) . (ICT . 1121), dlSl), gH), g31), gH) into (IXTOl) and arrive at the 
formula ()2.5|) for Cg^'^'. The proof is complete. 

5. Examples 

We shall now apply our results to obtain the expansion of the first eigenvalue for 
different domains, and compare the values obtained with a numerical approxima- 
tion. We are indebted to Pedro Antunes for carrying out the necessary numerical 
computations. 

We have chosen five examples illustrating several possibilities for the functions 
H, h+ and The first three correspond to algebraic curves, namely the circle, 
the lemniscate and the bean curve. For these h+ = h-, and the fit between the 
four-term asymptotic expansion and the numerical approximation is very good: for 
the lemniscate and the bean curve the error is always below 2% and for the disk 
the maximum error is around 5%. 

The last two examples have 7^ /i-, and the second one is nonconvex. Here we 
see that the error becomes much larger - in the last example it can go up to 50%. 

Example 1 (Disk). Consider the disk centred at (1/2,0) and radius 1/2 for which 
we have H[xi) = 2/i_|_(a;i) = 2h-{xi) = 2{x — x^)^!"^ . The maximum of H occurs 
at x\ = 1/2 and we obtain the expansion given by p.ip . Comparing these results 
with those of Daymond referred to in [j] we see that the error up to e equals one 
is maximal at one and is around five per cent. 

Example 2 (Lemniscate). Consider the lemniscate defined by 

( 2 . 2\2 _ 2 2 
\Xi + X2) ~ Xi ~ X2- 

In this case we have 

1/2 

H{xi) = 2h+{xi) = 2h^{xi) = 2 



l-xl + l{l + 8xlf' 



and the maximum of H is now situated at ^3/(2^2). This yields 
2\/37r 97 / 593 VStt' 



27r2 



Ai(e) = ^ + + — + \—^ + —-]e + Oie'), as e ^ 0, 



e 



24 



I 64%/37r 
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Example 3 (Bean curve) . As a third example we consider the quartic curve defined 
by xi{xi — l){xl + X2) + X2 ~ 0- We now have 



H{xi) = 2h+{xi) = 2h^{xi) = (2x1)^/^ 1 - xi + (1 - xi)^^^ (1 + 3x1)^/^ 

The maximum of H is situated at xi — 2/3 and we obtain 

, . , 97r2 3VT5TT 127 / 24229 bVbA 2, 

Ai e = n+— ^ + ^+\ ^ + ^K + C'e ' as e ^ 0, 

16e2 8e 40 \^600\/T57r 16^3/ 

Example 4 (Convex, with 7^ Let /i+(a;i) = sin(7ra;i) and h^{xi) = 7r(l — 

a;i)/2, yielding iJ(a;i) = sin(7ra;i) + 7r(l — a;i)/2, which has a maximum at = 1/3. 
The expression for the eigenvalue asymptotics now becomes 



1/2 



(3\/3 + 27r)2e2 (3\/3 + 27r)3/2 



9(27 + 6737r + IGtt^) _ 19 ' 
16(3^3 + 27r)2 ^ 

4^ 132 73 , 180 77r , 54657r^ , J72577r^\ 3~^/V^e 

3^/22304 31/^11664 ^ (3v^+ 2^)^/2 

+0(6^), as e ^ 0. 

Example 5 (Non convex). Let ft,+ (xi) ~ l + sin(77ra;i/2) and h-{xi) = 77r(l — a;i)/4, 
yielding H{xi) = 1 + sin(77r.Ti/2) + 77r(l — xi)/4, which has its global maximum 
on (0, 1) at xi = 2/21. We have 

n 210941 1 79fi92 
Ai(£) = + i^^yoyz ^ ^ ^ 60.57G6e + 0{e^), as e ^ 0, 

e e 

where, for simplicity, we only presented the numerical values of the coefficients. 
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